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Abstract 

Although a direct interaction between WIMP CDM candidates and 
scalar Quintessence fields (sources of Dark Energy) poses severe phe- 
nomenological threats, yet it is possible that the presence of Quintessence 
produces profound deviations in the expansion rate of the early Uni- 
verse at the moment of WIMP decoupling, hence causing significant 
enhancements of the relic DM abundance. We consider the occurrence 
of such a phenomenon in a couple of physically interesting situations, 
kination and scalar-tensor theories of gravity where the scalar compo- 
nent plays the role of Quintessence. The consequences for supersym- 
metric DM candidates are briefiy discussed. 



1 Introduction 

It is known that in the last years of the past century a major surprise oc- 
curred in our knowledge of the Universe: the dominance of Dark Energy 
(DE) in the total energy budget of the Universe, together with a close value 
of DE and Dark Matter (DM) energy densities in the present epoch of the 
Universe evolution. Indeed, their ratio is roughly two today, while one may 
expect it to have taken values different from 0(1) even by orders of mag- 
nitude during the whole evolution of the Universe from the Big Bang until 
the present time. This "coincidence problem" may hint at a possible corre- 
lation between DE and DM. For instance, if one is to envisage DE as result- 
ing from an ultralight scalar field rolling down its potential (the so-called 
Quintessence explanation of DE), then one could naively think that such a 
correlation results from a direct interaction between the DM particles and 
the scalar Quintessence field. However, the extreme lightness of the scalar 
Quintessence field, whose mass should correspond roughly to the inverse of 
the Hubble parameter, i.e. O(10~^'^eV), and its time variation are at the 
basis of a host of severe phenomenological problems to be tackled (new long 
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range forces with violations of the weak equivalence principle, time varia- 
tion of the fundamental couplings, etc.). In this talk we wish to emphasize 
that the severe limitations we have to impose on such direct interactions 
between DM and DE do not prevent DE to have an indirect impact on DM. 
Namely, we wish to examine a couple of situations where the presence of 
quintessential DE in the early Universe largely influences today abundances 
of DM with significant changes on our prospects for direct and indirect DM 
searches. 

According to the standard paradigm, a particle species goes through 
two main regimes during the cosmological evolution. At early times it stays 
in thermal equilibrium, until the particle interaction rate F remains larger 
than the expansion rate H. Later on, the particles will be so diluted by the 
expansion of the universe that they will not interact anymore and H will 
overcome T. The epoch at which F = H is called 'freeze-out', and after that 
time the number of particles per comoving volume for any given species will 
remain constant. This is how cold dark matter particle relics (neutralinos, 
for example) are generated. 

As it can be easily understood, this scenario strongly depends on the 
evolution equation for H in the early universe, which is usually assumed to 
be radiation-dominated. However, as it was already noticed some time ago 

there is little or no evidence that before Big Bang Nuclesynthesis (BBN) 
it was necessarily so. Non-standard scenarios are then worth exploring. In 
particular, if we imagine that for some time in the past the Hubble parameter 
was larger than usually thought (for example, due to the presence of some 
other component, in addition to radiation), then the decoupling of particle 
species would be anticipated, resulting in a net enhancement of their relic 
abundance. 

A natural candidate for doing that is the Quintessence scalar, which 
is thought to consitute the Dark Energy fiuid dominating the present uni- 
verse. In most Quintessence models, the cosmological scalar is assumed to 
become the dominant component of the universe after a long period of sub- 
domination W, playing little or no role in the earliest epochs. However, this 
has not always to be the case, as we will show in the following. We will 
focus on the possibility of modifying the past evolution of the Hubble para- 
menter H, with the double aim of respecting all the post-BBN bounds for 
the expansion rate and of producing a mesurable enhancement of the dark 
matter particles relic abundance. In particular, we will report about two 
possible mechanism by which the past dynamics of the Quintessence scalar 
could significantly modify the standard evolution of the pre-BBN universe: 
an early "kination" phase and a "scalar-tensor" model. 
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2 Kination enhancement 



If we imagine to add a significant fraction of scalar energy density to the 
background radiation at some time in the cosmological history, this would 
produce a variation in H^, depending on the scalar equation of state tf^. 
li > Wr = 1/3, the scalar energy density would decay more rapidly 
than radiation, but temporarily increase the global expansion rate. This 
possibility was explicitly considered in Ref. where it was calculated that 
a huge enhancement of the relic abundance of neutralinos could be produced 
in this way. 

In a flat universe, a scalar field with potential V{(j)) obeys the equations 

4) + 3H^ + dV/d()) = ; = {d/af = 8ttp/3M^ . (1) 

For any given time during the cosmological evolution, the relative impor- 
tance of the scalar energy density w.r.t. to matter and radiation in the total 
energy density p = pm + Pr + depends on the initial conditions, and is 
constrained by the available cosmological data on the expansion rate and 
large scale structure. If the potential V{(j)) is of the runaway type, the initial 
stage of the scalar evolution is typically characterized by a period of so-called 
'kination' 2 during which the scalar energy density p^ = <j)'^/2 + V{(p) is 
dominated by the kinetic contribution = (p^ /2 ^ V {(/)), giving = 1. 
After this initial phase, the field comes to a stop and remains nearly con- 
stant for some time ('freezing' phase), until it eventually reaches an attractor 
solution 0. 

Then, if we modify the standard picture according to which only radia- 
tion plays a role in the post-inflationary era and suppose that at some time t 
the scalar contribution was small but non negligible w.r.t. radiation, then at 
that time the expansion rate H{t) should be correspondingly modified. Dur- 
ing the kination phase the scalar to radiation energy density ratio evolves 
like Pifi/pr ~ fl-3(?i>0-u)r) _ g^j^^ gQ ^Yie scalar contribution would rapidly 
fall off and leave room to radiation. In this way, we can respect the BBN 
bounds and at the same time keep a significant scalar contribution to the 
total energy density just few red-shifts before. The increase in the expansion 
rate H due to the additional scalar contribution would anticipate the de- 
coupling of particle species and result in a net increase of the corresponding 
relic densities. As shown in [2], a scalar to radiation energy density ratio 
P(j)/ Pr — 0.01 at BBN would give an enhancement of the neutralino codensity 
of roughly three orders of magnitude. 

The enhancement of the relic density of neutralinos requires that at some 
early time the scalar energy density was dominating the Universe. This fact 
raises a problem if we want to identify the scalar contribution responsible 
for this phenomenon with the Quintessence field jl]. Indeed, the initial 
conditions must be such that the scalar energy density is sub-dominant at 
the beginning, if we want the Quintessence field to reach the cosmological 
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attractor in time to be responsible for the presently observed acceleration 
of the expansion [2]. For initial conditions p^f) ^ pr we obtain instead an 
'overshooting' behavior: the scalar field rapidly rolls down the potential and 
after the kination stage remains frozen at an energy density much smaller 
than the critical one. However, as shown in 5 , more complicated dynamics 
are possible if we relax the hypothesis of considering a single uncoupled 
scalar. The presence of several scalars and/or of a small coupling with 
the dark matter fields could modify the dynamics in such a way that the 
attractor is reached in time even if we started in the overshooting region. 

Consider a potential of the form V{(i)i,(j)2) = M"+*^ {(pi(p2)~'^^'^ , with M 
a constant of dimension mass. In this case, the two fields' dynamics enlarges 
the range of possible initial conditions for obtaining a quintessential behavior 
today. This is due to the fact that the presence of more scalars allows to 
play with the initial conditions in the fields' values, while maintaining the 
total initial scalar energy density fixed. Doing so, it is possible to obtain 
a situation in which for a fixed in the overshooting region, if we keep 
initially 4>i = 4>2 we actually produce an overshooting behavior, while if we 
choose to start with (pi ^ (j)2 (and the same p^") it is possible to reach the 
attractor in time. 

Suppose, instead, that the Quintessence scalar is not completely decou- 
pled from the rest of the Universe. Among the possible interactions, two 
interesting cases are the following: 

Vi, = bH^<i? or Vc = cp^4> (2) 

If we add Vh oi Vc to V = M^^^^"", the potential will acquire a (time- 
dependent) minimum and the scalar field will be prevented from running 
freely to infinity. In this way, the long freezing phase that characterizes 
the evolution of a scalar field with initial conditions in the overshooting 
region can be avoided. A more detailed discussion, together with numerical 
examples, can be found in Ref. j^. 

3 Scalar-tensor enhancement 

A different possibility arises if we consider Quintessence models in the frame- 
work of scalar-tensor (ST) theories of gravity (see and references therein). 
These theories represent a natural framework in which massless scalars may 
appear in the gravitational sector of the theory without being phenomeno- 
logically dangerous, since they assume a metric coupling of matter with 
the scalar field, thus ensuring the equivalence principle and the constancy 
of all non-gravitational coupling constants _6 . Moreover a large class of 
these models exhibit an attractor mechanism towards GR 0, that is, the 
expansion of the Universe during the matter dominated era tends to drive 
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the scalar fields toward a state where the theory becomes indistinguishable 
from GR. 

ST theories of gravity are defined, in the so-called 'Jordan' frame, by 
the action 

Sa = -^j d^xy/^ [^""R + 4a;($)5^'^a^$9,$ - 4F($)] . (3) 

The matter fields are coupled only to the metric tensor g^^, and not to 
i.e. Sm = Sm[^m, g^j.u]- Each ST model is identified by the two functions 
a;($) and V{^). The matter energy- momentum tensor is conserved, masses 
and non-gravitational couplings are time independent, and in a locally in- 
ertial frame non gravitational physics laws take their usual form. Thus, the 
'Jordan' frame variables and $ are also denoted as the 'physical' ones 
in the literature. By means of a conformal transformation, 

g^,=A\y,)g^, , = 8ttM^ A'^ip) (4) 

with 

a2(^) = dlogA{ip)/dip = l/(4a;($) + 6) , (5) 

it is possible to go the 'Einstein' frame in which the gravitational action 
takes the standard form, while matter couples to ip only through a purely 
metric coupling, 

Sm = Sm[^m,A\^)g^,] . (6) 

In this frame masses and non-gravitational coupling constants are field- 
dependent, and the energy-momentum tensor of matter fields is not con- 
served separately, but only when summed with the scalar field one. On the 
other hand, the Einstein frame Planck mass M^, is time-independent and 
the field equations have the simple form 

R^. - yf..R = t;^Jm^ + v/m2 , M^d^ip + dv/dip = -a{p^)T/V2 . 

(7) 

When 0(99) = the scalar field is decoupled from ordinary matter and the 
ST theory is indistinguishable from ordinary GR. The effect of the early 
presence of a scalar field on the physical processes will come through the 
Jordan- frame Hubble parameter H = dlogd/df: 

H = H{l + a{^)p')/A{^), (8) 

where H = dloga/dr is the Einstein frame Hubble parameter. A very at- 
tractive class of models is that in which the function a{<p) has a zero with a 
positive slope, since this point, corresponding to GR, is an attractive fixed 
point for the field equation of motion 8 . It was emphasized in Ref. |H] 
that the fixed point starts to be effective around matter-radiation equiva- 
lence, and that it governs the field evolution until recent epochs, when the 
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Quintessence potential becomes dominant. If the latter has a run-away be- 
havior, the same should be true for a{ip), so that the late-time behavior 
converges to GR. For this reason, we will consider the following choice , 

A{ip) = 1 + Be-^'f , a{^) = -/?Se-^^/(l + Be'^"^) , (9) 

which has a run-away behavior with positive slope. 

In Ref. it was calculated the effect of ST on the Jordan-frame Hubble 
parameter H at the time of WIMP decoupling, imposing on the parameters 
B and /? the constraints coming from GR test, CMB observations and BBN. 
Computing the ratio H/Hq^ at the decoupling time of a typical WIMP 
of mass m = 200 GeV, it was found that it is possible to produce an en- 
hancement of the expansion rate up to O(IO^). As a further step, it was 
performed the calculation of the relic abundance of a DM WIMP with mass 
m and annihilation cross-section (cJann^)- The effect of the modified ST 
gravity enters the computation of particle physics processes (like the WIMP 
relic abundance) through the "physical" expansion rate H defined in Eq. Q . 
We have therefore implemented the standard Boltzmann equation with the 
modified physical Hubble parameter H: 

dY/dx = -s{a,n^v){Y^-Y^l)/H X (10) 

where x = m/T, s = (27r^/45) hi,{T) is the entropy density and Y = nj s 
is the WIMP density per comoving volume. 

A numerical solution of the Boltzmann equation Eq. (jlOj) is shown in 
Fig. n for a toy-model of a DM WIMP of mass m = 50 GeV and constant 
annihilation cross-section (ciann^^) = 1 x 10^^ GeV^^. The temperature 
evolution of the WIMP abundance Y{x) clearly shows that freeze-out is 
anticipated, since the expansion rate of the Universe is largely enhanced by 
the presence of the scalar field ip. This effect is expected. However, we note 
that a peculiar effect emerges: when the ST theory approached GR (a fact 
which is parametrized by ^(y?) ^ 1 at a temperature T,^, which in our model 
is 0.1 GeV), H rapidly drops below the interaction rate F establishing a short 
period during which the already frozen WIMPs are still abundant enough 
to start a sizeable re-annihilation. This post-freeze-out "re-annihilation 
phase" has the effect of reducing the WIMP abundance, which nevertheless 
remains much larger than in the standard case (for further discussion on 
this aspect see (7]). The amount of increase in the relic abundance which 
is present in ST cosmology is shown in Fig. |2l The solid curve refers to an 
annihilation cross section constant in temperature, i.e. (cxann^^) = Oi while 
the dashed line stands for an annihilation cross section which evolves with 
temperature as: {(JwrnV) = hjx. In the case of s-wave annihilation the 
increase in relic abundance ranges from a factor of 10 up to a factor of 400. 
For a pure hjx dependence, the enhancement can be as large as 3 orders of 
magnitude. Needless to say, such potentially (very) large deviations entail 
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B, Catena. N. Fornengo. A Masiero, M. Pietroni, F. Rosali (2004) 




x = m /T 



Figure 1: Numerical solution of the Boltzmann equation Eq. (jlUj) in a ST 
cosmology for a toy-model of a DM WIMP of mass m = 50 GeV and con- 
stant annihilation cross-section (dann^^) = lxlO~^ GeV~^. The temperature 
evolution of the WIMP abundance Y{x) clearly shows that freeze-out is an- 
ticipated, since the expansion rate of the Universe is largely enhanced by 
the presence of the scalar field ip. At a value x = m/T^ a re-annihilation 
phase occurs and Y{x) drops to the present day value. 



new prospects on the WIMP characterization both for the choice of the 
CDM candidates and for their direct and indirect detection probes [7| . 
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